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^ ! Abstract. 

Q"^ I We describe applications of (perturbed) conformal field theories to two-dimensional disor- 
dered systems. We present various methods of study : 

(i) A direct method in which we compute the explicit disorder dependence of the corre- 
al lation functions for any sample of the disorder. This method seems to be specific to two 
r-| ! dimensions. The examples we use are disordered versions of the Abelian and non-Abelian 
WZW models. We show that the disordered WZW model over the Lie group Q at level k is 
K/i ; equivalent at large impurity density to the product of the WZW model over the coset space 
' Q'" /G at level {—2h'") times an arbitrary number of copies of the original WZW model. 



(ii) The supersymmetric method is introduced using the random bond Ising model and 
the random Dirac theory as examples. In particular, we show that the relevent algebra is 
the affine OSp{2N\2N) Lie superalgebra, an algebra with zero super dimension. 

(iii) The replica method is introduced using the random phase sine-Gordon model as 
example. We describe particularities of its renormalization group flow. 

(iv) A variationnal approach is also presented using the random phase sine-Gordon model 
as example. 



"'^ Lectures presented at the '95 Cargese Summer School on "Low dimensional application of quantum field 
theory" . 

^Member of the CNRS 

^ Laboratoire de la Direction des Sciences de la Matiere du Commisariat d I'Energie Atomique. 
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1 Introduction. 



These lectures will probably look too elementary to experts in disordered systems. It will 
perhaps be more appropriate as an introduction to this field for conformal field theorists. 
The plan of these notes is the following: 

1- Introduction. We introduced the types of problems we will be concerned with. We 
mention a few open problems; the most challenging being the description of the quantum 
Hall transition. 

2- A few basic techniques. In this section we have gathered the basic tools which are 
usually used. This includes the supersymmetric method and the replica method. 

3- The direct method. There, we show that due to specific properties of two-dimensional 
field theories, it is (sometimes) possible to compute directly the disorder dependence of the 
correlation functions and then to average over the disorder. Examples are provided by WZW 
models coupled to random vector potentials, a particular case being the massless Dirac model 
coupled to a random vector potential. We discuss the Abelian and non-Abelian cases and 
show in both cases that the disorder can be factorized. At large impurity density, the dis- 
ordered WZW model is conformally invariant being equivalent to the product of the WZW 
model over the non-compact space Q'" /Q times an arbitrary number of copies of the pure 
version of the original WZW model. 

4- The supersymmetric method. We take the random Ising model as the simplest 
example for illustrating the supersymmetric method. We show that it can be formulated 
as a Gross- Neveu type model but over the Lie superalgebra OSp{2N\2N), and we use this 
formulation to recover its properties usually obtained by the replica method. We discuss its 
connection with the random Dirac theory: a model which has been introduced for describing 
the quantum Hall transition. These examples illustrate the fact that non-unitrary conformal 
field theories based on affine Lie superalgebras with zero superdimension provide candidates 
for critical theories for gaussian disordered systems. 

5- The replica method. We take the random phase sine-Gordon model as example. 
Its large distance behavior in the low temperature phase is still controversial. We present 
the approach based on the renormalization group and the symmetric replica trick. We also 
discuss a large version of this model. 

6- The variational method. We again use the random phase sine-Gordon model to 
illustrate this method. This approach consists in finding a good variational ansatz at fixed 
disorder and then averaging. It provides exact lower bound to the free energy. 

Besides stochastic differential equations, cf. Polyakov's lecture in this volume, there are 
at least two classes of disordered problems which we can study. 

i) The first ones concern statistical systems with random coupling constants, cf eg. 
Imagine having a physical system with two kinds of degree of freedom with two very differ- 
ent relaxation times. The degrees of freedom, called spin variables, which have the shortest 
relaxation time will thermalyze well before the other degrees of freedom, which we call im- 
purities. The impurities are also called quenched variables because they are not in thermal 
equilibrium with the spin variables. We consider a large number of such physical systems 
with different realizations of impurity configurations. Since the spin variables are in ther- 
mal equilibrium but not the impurities we are interested in the average over the impurity 
configurations of the free energy of the thermalyzed spin systems. 

A classical example is provided by the Ising model with random bond coupling constants. 
The thermalyzed spin variables (jj, defined on the vertices of a lattice, take two values ai = ±. 
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The effect of the impurities is represented by the randomness of the interaction Jij which 
couple neighbourhood spins. At fixed disorder, its partition function is : 

Z[J] = ^ exp - ^ Jijaiaj 

A probabihty distribution P[J] is given to the random couphng constants. The average ther- 
mal properties will then be given by the average of the logarithm of the partition function : 

In the continuum limit and close to the critical point, ie. in its scaling limit, it is known 
from prehistory that the Ising model is described by a massive real Majorana fermion with 
mass m ~ r = with Tq the critical temperature. In presence of randomness in the 

bond interaction, the mass becomes a random function of the space position. Therefore, at 
fixed disorder, the action describing the random bond Ising model in its continuous limit is : 

— [^i0'^ + m(x)\I^\I'j 

with a Majorana fermion. The random mass m(x) is coupled to the energy operator as 
it should be in order to represent the randomness of the bond interaction. This action is of 
the form : 

S[g{x)] =S. + J (fx g{x)^{x) (1) 

where 5*,, is the action at the renormalization group fixed point describing the pure system, 
and $(2;) a scalar field. As formulated in eq.(|l]), the problem can also be interpreted as a 
random field problem. 

The Ising model with random bond interactions provides one of the simplest example 
of disordered system in 2d 0, ^, ||. We will study it in details in the following (using 
the supersymmetric method). Besides the disordered WZW models, we will also study 
another example (using the replica method) : the random phase sine-Gordon model. It 
has been introduced for describing a large variety of two dimensional disordered physical 
systems H, 0, §] . Contrary to the random bond Ising model which is now well understood, 
there is up to now no concensus on the large distance behavior of the random phase sine- 
Gordon model. We will present the various suggestion which have been made concerning 
this behavior. 

ii) The second class of problems deals with particles moving in random potential, cf eg. 
[0]. The random potential then represents the impurities that the particles, e.g. the electrons, 
encountered while travelling in a metal. The property of the medium will be characterized by 
the behavior of the electrons in this disordered surrounding. For example, its conductivity 
properties will depends on whether the averaged wave functions are localized or not. These 
problems are generically called localization problems. 

We are thus interested in the averaged properties of hamiltonians of the form : 

H = Ho + V 

where Hq is a fixed hamiltonian, and V a random perturbation. A standard example is : Hq 
a pure kinetic term, Hq = —d^ and V = V{x) a potential representing the interaction with 
the underlying medium. 
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Quantities of interest would be e.g. the averaged density of states which is given by the 
average of the trace of the Green function : 
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More generally, we are interested in computing averages of products of Green functions. 
Indeed, computing the conductivity tensor a^y using the Kubo formula requires computing 
the average of the product of an advanced and a retarded Green function. 

Let me introduce a recent localization/delocalization problem (which actually was the 
motivation for these lectures). It concerns the transition between the plateaux of the integer 
quantum Hall effect, cf |P, |lOl and references therein. Without going into any details, we 
recall that the quantum Hall effect is characterized by the quantization of the transverse 
Hall resistance Rxy to inverse integer values in unit of (h/e^). This integer quantization 
remains valid when the magnetic field varies on some domain, which are called plateaux. 
On these plateaux the longitudinal resistance Rxx vanishes, but it becomes non-zero in the 
regime separating two plateaux. The values of the resistance R^y and R^x as a function of 
the magnetic field is represented in Fig.l. 




Fig.l : The transverse and longitudinal Hall resistances Rxy and Rx 
as a function of the magnetic field B. 



The integer quantum Hall effect (IQHE) is correctly described by considering a 2d gas 
of non-interacting electrons subject to a transverse magnetic field. As is well known, the 
energy spectrum is then given by the Landau levels, which are highly degenerate. In a pure 
system without any impurities, the Hall resistance Rxy will be a step-function as a function 
of the Fermi level. The plateaux occur when the Fermi level of the electrons gas is between 
two Landau levels. 

However, since the degeneracy of a landau level is field dependent, (this degeneracy is 
0/00 where is the magnetic flux through the system and 0o = h/e the flux quantum), the 
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Hall resistance in a pure system is a linear function of the magnetic field B. More precisely 



j^pure ( —— 



h . 
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where Ne is the number of electrons in the system and z/= ^ is called the filling factor. 
The explanation of the plateaux for R^y as a function of the magnetic field B and not of the 
Fermi level requires considering the effect of the impurities |1TT| , [12|. The impurity potential 



splits the degeneracy of the Landau levels, which now form energy bands. The states whose 
energy is in the border of these bands are localized while those which are in the center of the 
bands are delocalized. The localized states then serve as reservoir for the electrons which 
stabilized the Fermi level between two Landau bands while the magnetic field varies. This 
stabililization of the Fermi level give rise to the plateaux. It is expected that the delocalized 
states are present only at the center of the band; that is there is not a band of delocalized 
states. The presence of delocalized states is necessary for having a non-zero conductivity. 




Fig. 2 : The density of states in two consecutive Landau bands. 



Between two plateaux the Fermi level is (almost) in the middle of a Landau band. Since 
the delocalized states are in the middle of this band, the transition between two plateaux of 
the integer quantum Hall effect is thus a transition between localized and delocalized states. 
(De)localized states are charactized by a (in)finite localization length. Near this transition, 
this localization length C,{E) behaves as follows : 

^{E) \E - E,]-'' (3) 

with Ec the band center energy. The density of states is non-singular at the transition point 



1^. Various mesurements, as well as numerical simulations, gave access to the exposant 
/ : z/ ~ 2.3 ± 0.1 ~ 7/3 0. The behavior (H) has to be compared to the standard 



Anderson localization in two dimensions which predicts that any tiny amount of disorder 
will localize all the states. 

What is the field theory describing the quantum Hall transition is still an open question, 
cf ref. |TD| for a recent discussion. The problem consists in determining the appropriate in- 



frared fixed point, and therefore seems to be a "simple" exercise in conformal field theory. 
Unfortunalely, this infrared fixed point is a strong coupling problem in all the models intro- 
duced so far; it is therefore difficult to attack. A sigma model approach was introduced in 



5 



ref.[16, 17|. A very elegant network model was developped in ref.[lS]. A supersymmetric spin 



chain, which arises as an anisotropic limit of the network model, was introduced in ref.jlO 



A model closely related to the random bond Ising model, namely the random Dirac model, 
has also been proposed for describing the IQHE localization/delocalization transition 
We will present it in section 4. 
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2 A few basic techniques. 

These problems can be analyzed with similar techniques. The two standard ways of studying 



them use either the replica trick, cf eg or the supersymmetic method [20|. In this section 
we present a very brief introduction to these methods. (The similarity of the techniques used 
to study disordered statistical model and localization problems is apparent at a first naive 
level, but more refined studies of both subjects reveal their differences.) 

As it will become clear below, the supersymmetric method applies to system which are 
gaussian at fixed disorder, while the replica trick has a larger domain of applicability since 
it is not restricted to free theory. Perturbatively the replica and supersymmetic methods 
are equivalent. However, a few non-perturbative results have been obtained using the su- 
persymmetric method. Clearly, once these disordered problems have been reformulated in 
a field theoretical languages, all the standard field theoretical methods (e.g renormalization 
group, mean field approximation, variational approximation, large N developments, etc..) 
are avalable. 

In the following sections we will illustrate these methods on two examples : the random 
bond Ising model and the random phase sine-Gordon model. 

•) The supersymmetric method. 

The supersymmetric method |^ applies to models which are gaussian at fixed disorder. 



We will present it in the random Hamiltonian problem. As above let us consider a random 
hermitian hamiltonian 

H = Ho + V, 

where Hq is a fixed hamiltonian and V a random potential. We are interested in computing 
averages of product of advanced or retarded Green functions : 



He — iej \ He' + 

with e ^ 0^. We have set He = H — E. The supersymmetric method is based on a gaussian 
integral representation of the Green functions. E.g.: 



1 1 



He - ie Det[HE - ie] 



dtp tpip* exp {±iij*{HE - ie)ip) (4) 



with ip complex grassmanian variables. In eq.@, we have the freedom to choose the ± sign. 
There are no convergence problem in this fermionic integral. To be able to average over 
the disorder, we need to re-exponentiate the determinant Det[HE — ie]. The inverse of this 
determinant can be written as a gaussian integral over complex bosonic variables : 



Det[HE - ie] 



exp{—i(j)*{HE — ie)(j)) (5) 
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The prefactor (— z) in the exponential is choosen by requiring the convergence of the gaussian 
integral. Gathering these gaussian integrals, we obtain the supersymmetric representation 
of the Green functions : 

J dipdcp ipi;* exp {-ii;*{HE - ie)ip - i<p*{HE - ie)0) (6) 
j dipd(f) 'ipn'ilj*j^'iljAi'*Aexp {-iip}^{HE - ie)ipR - - «e)0ij 

In the equation representing product of two Green functions, the indices A, R refers to the 
advanced or retarded Green functions. Note that due to the convergence requirement, the 
sign in front of the (ze) factor differs in the advanced and retarded sector. More generally, 
products of N Green functions are represented by integrals over N couples of superpartners 
(-0, (f) with appropriate choice of the sign factor according to the convergence condition. 

In the form (^, the average over the disorder can be easily done. Assume for simplicity 
that the random potential has a gaussian distribution : 

P[V] = exp (^-^tr{V^)^dV 

Consider for example the average of one Green function. The term coupled to the random 
potential is ^trlV"^) + itr{y{(j)(j)* — ipip*)). After integration over V , the averaged Green 
functions becomes : 

{jf—^ = / di'd(p iIjV exp {-Seff) (7) 

with 

Seff = triHo -E- it)^ + 20* (/7o -E- 2e)0 - |tr (00* - ^^f 

The bosonic and fermionic sectors are coupled by the disorder. This effective action is easily 
generalized when considering products of higher number of Green functions. 

Notice that if if is a quantum hamiltonian for a physical system in d spacial dimensions, 
ie. a ((i + 1) system, than the supersymmetric action describes a d dimensional quantum 
field theory. In particular, since the quantum Hall effect is a (2 + 1) system its localiza- 
tion/delocalization transition would be described by a 2d field theory. 

•) The replica method. 

Let us now recall briefly the starting point of the replica method, cf We are interested 
in computing the average of the free energy, or of correlation functions. Since the free energy 
is proportional to the logarithm of the partition function Z[J\, we have to compute the 
average of logZ[J]. The replica method is based on the identity : 

log Z\J] = lim ^"^'^^ ~ ^ = -^Z"[J] (8) 

In other words, we compute the quenched averages of one disordered system by introducing 
n copies of this system, in all of which the random variables take the same values, and then 
computing the averaged partition function of this new system at = 0. 



He - ie 



He — it) \He' + ie 



...etc. 
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This method leads immediatly to a criteria selecting the systems in which the disorder is 
relevent. This criteria is called the Harris criteria. 

Let us describe a disordered system in the continuum limit by the action : 

S[g{x)] =S, + Y.(d?x g''{x)^a{x) (9) 

where is the action at the renormalization group fixed point describing the pure system 
at its critical point, and ^a{x) are scalar fields. When the coupling constant g°'{x) are 
independent of x and not random, the action describes the behavior of the statistical 
model near the critical point in its scaling limit. In this case the field ^a{x) is relevent, and 
therefore will infiuence the large distance behavior of the pure system, if its scaling dimension 
is less than two : dim($a) < 2. This is changed when disorder is present. 

Assume for simplicity that the random variables g"'{x) are gaussian variables with one 
and two-point functions : 

gHx) = and g''{x)g^{y) = (5"Va5(x - y) (10) 

Before averaging over the disorder, the partition function is Z[g{x)] = J DX exp{—S[g{x)]) 
and the n replicated one becomes : 

Z^'igix)] = fY[DXr exp f- E - E / d'x g\x) E K{x)\ 

r \ T a r / 

where the index r runs from 1 to n, SI is the action for the r^^ replic and $a(a;) is the field 
^a{x) in the r*^ copy. To have introduced the replicated copies allows us to easily average 
over the disorder. Assuming the gaussian measure (|TUp for g"'{x), we obtain : 

= [UDXr exp i-S^ff) (11) 

with 

Seff = E ^: + ^ E / rf^x E Kixmx) (12) 

r ^ a r,s 

The replica are now coupled by the disorder. The field perturbing the uncorrelated action 
J2r is quadratic in terms of the $a(x)'s. It will be relevent if 2 dim ($a(a;)) < 2. (A more 
careful study is needed to analyse the short distance singularity in $^(a;)$^(x) when r = s). 
Hence : 

dim^a{x) < 1 ^ relevent disorder. 

This is the Harris criteria in two dimensions. Note that for a field to be relevent in a quenched 
system its dimension has to be half what is required in the pure system. If these dimensions 
depend on parameters the Harris criteria provides a first glance to the phase diagram. 
The quenched correlation functions are defined by : 
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The replica method can also be used to compute the averaged correlation functions. It 
provides a way to re-exponentiate the inverse partition function. The rules are the following : 



{0{x)0{y)) = (0^(x)0^(y))U=o (13) 



{0{x)){0{y)) = {O'^ix)0^iy))\n=o, r^s 
... etc .... 

In particular, the average of the connected Green functions are represented by the correlation 
function of - X^r 0'^{x). The simplest way to prove eq.([T3|) consists in introducing sources in 
the quenched partition function. 

3 The direct method : the random vector potential 
model. 



As a warm up we introduce a very simple gaussian model, which evidently can be completely 
solved. (We will later use a few of its properities in our study of the more interesting random 
phase sine-Gordon model.) 

In two dimensions it is sometimes possible to compute exactly field theory partition 
functions and correlation functions with sources. In this section, we also described how this 
property can be applied to very particular disordered 2d systems. The example we choose 
consists in a non-Abelian generalization of the gaussian model based on the WZW models. 
It provides one of the rare disordered models which can be solved directly without relying 
on the supersymmetric or replica method. This solution relies on the fact that the disorder 
dependence can be factorized. 

•) The gaussian model. 

The bosonic form of its action is defined to be : 



cPx (\ 



TX 



(14) 



where A^ is a random quenched variable with measure : 







<" d?x 






<" (fx 


exp 
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= exp 




AzA- 



(15) 



In eq.(|T^ we introduced the complex coordinate z = x + iy, z = x — iy on the plane. 

Since it is a gaussian model, it can be easily solved without using the replica trick or 
the supersymmetric formulation, but by directly computing the quenched averages. Let us 



introduce the Hodge decomposition of A^ 
in the measure (p!5| ) : 



P[A] = P[^;ri] = exp 



A^ = d^^ + e^uduV- The fields ^ and i] decouple 



(16) 



'-±1 


' (Px / 


2a J 


47 V 
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The action 5*0 [$, A] is independent of ^ and therefore the field ^ is irrelevent. This fact was 
expected since the field ^ represents a pure gauge whereas the physically relevent quantity 
is the field strengh F = e^yd^A^ = {d^d^)rj. Moreover, the field r] can be absorbed into a 
translation of $ : 

S,[^,A, = e,,d.r^] = S,{^ + ir^,A = Q]+ J ^^{d,v)'. (17) 

The fact that the field rj can be absorbed into a shift of $ does not mean that the 
quenched correlation functions are identical to those in the pure system. Using eq.(|17D, we 
have 

n n 

It factorizes into the product of the correlation functions in the pure system times a simple 
function of the impurities. However, the average of this function is not irrelevent since the 
variables rj have long-range correlations : ri{x)ri{y) = —na{dzdz)~^{x,y) = — aloglo: — y\'^. 
In particular, it changes the values of the critical exponents. 

Averages of products of correlation functions can then be computed . One has : 



n m 

n m 

where ((■■■)) refers to the ?7-correlation functions with the free field measure (|T^). All the 
correlators are therefore given by gaussian integral. We get : 



n m n,m 

X n K - ^nf • n IVm - ym'^"^^-' 

n<n' m<m' 

Clearly, conformal invariance is unbroken in the random abelian case. The dimensions A of 
the vertex operators exp(m$(x)) in the quenched theory are : 

AO A a nnP' 

quenched pure ^ ^ " 

It is interesting to look at the connected Green functions of the current = id<^. It turns out 
that all these quenched connected correlations are chiral; in particular, they are holomorphic 
and only depend oti z = x + iy. For example, we have : 



{id^{z)) = 0, 



{id^z)id^w))conn. = 2TT[d=. 

o / z,w yz w)^ 
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More generally, one verifies that : 



{^^^Z,)■■■^^^Zp)),onn. = {td^z^)---td^zp))A=o- (19) 

In other words, the averages of these connected correlation functions are unaffected by the 
disorder. This properity is not true for the average of disconnected correlation functions of 
J, = id^. 
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A more complete derivation of these Ward identities which follows from a U{1) symmetry 
is presented in the section devoted to the study of the random phase sine-Gordon model. 
It is an interesting exercise to decipher the operator product algebra of these models. 

•) The non-Abelian case : factorization. 

We now turn to the non-Abelian case which provides an exceptional quenched model 
which can be solved directly, without using any replica or supersymmetric method. This is 
done by factorizing the disorder dependence. We derive it directly from the definition of the 
quenched correlation functions. This result was also found in ref.|]21| using the supersym- 
metric method, (it appears while we were preparing these notes). 

The non-Abelian model appears naturally when generalizing the previous model, after 
fermionization, to a iV-component model. Indeed, using standard bosonization/fermionization 
rules, \1/7^\1/ = e^^^d^^, and \l/i^\E' = we can rewrite the action (p^if ) in an equivalent 

fermionic form : 



^0 = J^{n0+m 



It describes a Dirac fermion coupled to a random magnetic field. 

The non-Abelian generalization describes N massless Dirac fermions minimally coupled 
to a random non-abelian gauge field P2| , |2T|. Let us introduce the components of the 
fermions {4'+, '^+) and '^-■,3)1 j = '^i ' ' ' 1 N . The action is defined as : 



(20) 



where A^ = iJ2aA'^{t°')i, with (A")* = A^, is the gauge field. Here the hermitian ma- 
trices form the A^-dimensional representation of SU{N). We denote by f"''"^ the SU{N) 
structure constants : = if°'^H'^. The Dirac fermions take values in this A^- dimensional 

representation. The gauge field is assumed to be a quenched variable with the gaussian 
measure : 



P[A] = exp 



-- f—T.A:A^ 



(21) 



At fixed disorder the partition function is : 



Z[A] 



Dtpe 



-S[A] 



Detlia 



(22) 



It can be expressed in terms of the WZW action see below eq. 

The way to factorize the disorder consists in implementing a chiral gauge transformation. 
So we parametrized the vector potential by an element G of the complexified group SU (A^)*" 
as follows : 



A-=G-' {d,G) = -Al (23) 

This is always possible on the sphere provided the connexion A is regular enough. This 
parametrization is unique up to left multiplication by a constant group element. Notice 
that the number of degree of freedom is preserved by this transformation : in the A or G 
parametrization there are 2dimSU{N) degree of freedom. We then have ip-{d^ + A-)iIj+ = 
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{i>-G ^)dz{Gijj+) and ip_{dz + Az)i/j_^_ = {ilj_G*)dz{G ^ We can therefore absorbed the 

dependence on A by transforming the fermions as : 

^_ =%l)_G-^ ■ C = ^-G* 
$+ = G^P+ ; ^ = G-i*V^+ (24) 



This is the chiral gauge transformation. It maps the action (pOD into the free Dirac action for 
the gauge transformed fermions. We denote it SfreeW\- The Jacobian for the transformation 
ip ^ 'ip is non-trivial but given by the chiral anomaly : 



Dip Det[z0] 

The extra factor Det[i0] is irrelevent when computing the quenched correlation function but 
relevent in the evaluation of the conformal anomaly. The crucial point is that this Jacobian 
is equal to the partition function at fixed disorder. Therefore, the partition and the Jacobian 
simplifies when computing the correlation functions at fixed disorder, and we get : 

D-ip e-^^'--[^] ■ • • (25) 

This is the announced factorization property : the correlation functions at fixed disorder 
factorize into the product of a correlation function in the free Dirac theory times an explicitely 
known functional of the vector potential. 

We now have to average over the disorder. For this it is convenient, if not necessary, to 



change variable from A to G. The Jacobian is known, see eg. |25 



DA = DG exp ( 2h'' S^,^{GG*) ) (26) 

where DG is the Haar measure on SU{N)'^, S^zwiGG*) is the WZW action and is the 
dual Coxeter number, equals to is the SU{N) case. Since only the product GG* enters into 
the Jacobian, we decompose SU{N)'^ as {SU{N)'^ / SU{N)) x SU{N). For G this means : 

G = H U with H G SU{Nf /SU{N), U G SU{N) 

The Haar measure DG factorizes on the Haar measure on [SU {N)^ / SU {N)) and SU{N). 
Thus, we finally obtain : 

DA = DUDHexp (2/i'^5^,^(M*) ) (27) 

In eq. (pT]) the part related to the compact space SU{N) only involves the Haar measure, 
while the part related to the non-compact symmetric space {SU{N)'-^ / SU{N)) involves the 
WZW action. Although the WZW action seems to appears with a 'wrong' sign in eq. (|27|) , the 
non-compact sigma model is well defined since the metric on {SU{N)'-' / SU{N)) is negative 
definite. 



For an arbitrary value of the disorder strengh a we also have to include the factor (^TD 
into the measure for A. This factor couples U and H. But things simplify at cr = oo. For 
this value of the coupling constant, we have : 

P.=^[A] DA = DUDHexp {2h'' S^.^^HH*) ) (28) 
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The compact and non-compact sectors are now independent. Moreover, the field U of the 
compact sector completely decouples from the quenched correlation functions of the local 
spinless neutral operators as e.g. (V'-V'+)- Indeed, we have : 



■■■)a--- ,,=00 = j DHDi) e^h-s^.^{HH*)-Si,,AM jjjj* . . . ^29) 

Notice that a necessary condition for the cancelation of the ^/-dependence in eq.(^) is to 
consider SU{N) scalar operators. le. we have to sum over the color indices of the fermions 
and in order to multiply U and U* to get UU* = 1. 

The free Dirac action plus the non-compact sigma model are conformally invariant. This 
is clear for the free massless Dirac theory. The sigma-model on the non-compact space 
{SU{N)'^ / SU{N)) defined by the WZW action was studied in ref.[^ in connexion with 
path integral construction of the coset models. There it is was shown that it is a non-unitary 
conformal field theory which carries a representation of the affine Kac-Moody algebra of 
negative level k = —2h'". Therefore, the Virasoro central charge of the non-compact sector 
is : 

C-^-=^l^^§0^ = 2MraSUiN) = 2(N^-l) (30) 

The scaling dimension of the primary fields (j)R which belong to a representation R are given 
as usual by : 

^^"^ 2((-2/i'') + /i^) 2/1'' ^ ' 

with Cas{R) the Casimir operator in R. The negative value of this dimension is an echo of 
the non-unitary character of the theory. 

As an application of eq. (|29|) we can derive the scaling dimension of the operator 
in the quenched correlation functions. It is the sum of the scaling dimensions of the op- 
erator {ip-ip^) in the free Dirac theory plus the one of the operator (HH*) in the vector 
representation (□) in the non-compact WZW model. Explicitely : 



dim^=oo[{'ip-'ip+)] = 1-2 



Cas{n) _ ^ _ ^ 



2N) m 



This agrees with ref. [22, 21 



In summary, we have shown that at a = 00 we have the following equivalence, valid in 
the local spinless neutral sector : 



DiracA a=oo = {Free Dirac)^ x {SU {N)^ / SU {N)y, 



with M the number of correlations we are averaging. This equivalence is a direct consequence 
of the factorization property eq.(^,^) of the random vector potential model. 

•) Quenched current correlation functions. 

The previous factorization property is really useful only at infinite coupling a. Here we 
describe how averages of current correlation functions can be computed for any value of a. 
We also show that, as in the Abelian gaussian model, the averages of the connected current 
correlation functions are unaffected by the disorder [2^. 
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We are thus interested in computing the quenched average of the correlation functions 
(\l/7^t'^\l/). Exphcitely, their components are : 



of the currents Jf, 



(32) 



While these currents are conserved in the pure system, they are not conserved once the 
disorder has been turned on. However, as explained in ref. [^, the quenched theory still 
possesses a su{N) symmetry. The currents generating this symmetry are represented by the 
insertion of the following operators in the quenched averages : 



6 



TT- 



Although the fields J" are not conserved, their quenched correlation functions can be com- 



puted. This relies on the Polyakov-Wiegman (PW) formula [24| for the effective action W[A] 



defined in eg . (|2^) . It can be exactly computed by integrating its anomalous transformation 
under a chiral gauge transformation. Indeed, let = with = = 

^^^^^ be the generating functions of the connected current Green function in the pure 



TT 



system. They satisfy the anomalous Ward identities |24 



d,Gj+c^G, + [A,,G-] + [A-,G, 
d,G--djG, + [A,,G-]-[A-,G, 



0, 

2F,^[A], 



(33) 



with F^j[A] = O^Ay - djA^ + [A^, A^]. Eqs.(||) completely specify G^[A]. The solution of 
eqs.(p3D can written as : 



GM 
Gj[A] 



A. 



1 



dz + ad.Az 
1 



d.A^ 



Oz + ad.Az 



(34) 
(35) 



Here ad.Az denoted the adjoint action of Az. Notice that G'z[A] is local in Az but non-local 
in Az. 

Knowing explicitely the correlation functions for any impurity sample, it is a priori pos- 
sible to take the quenched average. Let us first concentrate on the average of the chiral 
correlation functions involving only currents of the same chirality; e.g. involving only J". 
Consider first the average of products of one-point functions. Since the quenched average is 
defined by A'^{x)A^{y) = aTr6'^^5^'^\x — y) and since {Jz)a are linear in Az, these quenched 
correlations can be computed using Wick's theorem applied on A. For example, the two-point 
and three-point functions are : 



1 



lab' 



ah 



{zi,Z2) (^1 - 
1 



3a^ 



(zi,Z2) 

+ {cyclic permutation) 

^ Jabc 

{Zi - Z2){Zi - Z3){Z2 - Z^) 



I) 

iz3,Z2). 



14 



More generally, the average of products of one-point functions is the sum of connected 
correlations which can be expressed in terms of the correlation functions of the pure system : 



{JH^i))a ■ ■ ■ (^m))aJ = Ma''~\j:-{z,) ■ ■ ■ Jr{zM))o (36) 

Here, (■ ■ ■)o denote the pure correlation functions. They are known exactly The relation 
between the quenched correlations and their connected parts is the usual one. 

More interesting is the average of products of correlations with insertion of the conserved 
currents since they encode the underlying symmetry algebra. We find : 

= {j:^{z,) . . . jr{zM)JT{y^i)JlKii) • ■ ■)o ■ ■ • {jr{wp)Ji-{ip) ■ ■ -y 

+ E('^r(^i) ■ ■ ■ J^\z,) ■ ■ ■ Jr{zM)JT{wi)JlKii) ■ ■ ■)o X (37) 

X {j:^{z,)JT{w,)Jt^{i,) ■■■Y--- {j^M-nn^p) ■ ■ -y 

+ ■■■ 

+ (Jr (^i)^'Hei) ■■■Y--- {JTi^i) ■ ■ ■ JT'{zM)Jr{wp)Jl-{ip) ■ ■ -Y 

The hatted fields have to be omitted. Here, we assumed that there is no insertion of one-point 
functions. The formula ( P?] ) is actually simpler in words : it is obtained by distributing the 
currents J'^^{zi), • • • , J^^^^zm) among the pure correlators in all possible way, each counted 
only once. Notice that all the chiral quenched correlation functions are purely algebraic, 
without any logarithmic correction. 

Eq.(^) with no insertion of JT" shows that the connected correlation function of the 
fields J" are unaffected by the disorder. 

From eq. (|37D we read the operator product expansion of the fields. The currents J'^ 
satisfy : 

J:\z,)J:-{z,) = + ^ J--iz2) + reg. (38) 

(^1-2:2)^ Z1-Z2 

Therefore, the quenched conserved currents satisfy the commutation relations of a Kac- 
Moody algebra, exactly as the currents in the pure system do. 

The operator product expansion between the conserved currents J'^ and the correlators 
{J^'iw,)J^^iw2)---)A are : 

rnoj - — - 

j:iz){j:^iw,)j:^iw2)---YA ~ _ ,2 ('^:'M---J?i^,)---)A 

+ E "^(^r (^1) ■ ■ ■ JA^^) ■■■rA + reg. (39) 

■ Z Wj 

These operator product expansions are unusual in conformal field theory with Kac-Moody 
symmetry. In particular, they imply that the fields (J^i(wi) 7^2(^2) ■ ■ -Ya ^lot associated 
to highest weight vector representations. We don't know to which category of representations 
they correspond to. 

Contrary to the chiral quenched correlation functions which are easy to compute, the 
averages of correlation functions involving fields of opposite chiralities are difficult, if not 
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impossible, to evaluate directly. This is due to the fact that the generating functions Gz[A] 
and ^^[A] are non-local in and Az respectively. A naive perturbative expansion is spoilled 
by untractable divergences. However, since the model is asymptotically free in the infrared. 



one can used renormalization group techniques to evaluate few two-point functions 1^3 



•) Generalization to any WZW model. 

The previous discussion generalizes easily to any WZW models. So let us consider a 



WZW model over the Lie group Q at level k. Its action is |3^, |26 



S^zM = k I -^tr{dzg~'d,g) + kT{g) 

J OTT 



where T{g) is the WZW topological term. We define the action of the disordered WZW 
model to be the WZW action in presence of sources. Namely : 



S^zwig; A) = S^zwig) - ^ / ^ {tr{AzJ- + A-Jz) + tr{gAzg-^Az) - tr^AzA-) ) 

(40) 



where Jz, Jz are the WZW currents and Az, Az the random variables. We assume that the 
probability distribution P[A] is given by eq. fplf ). 

As a functional of g and A the action (HO) satisfies the "Polyakov-Wiegman" relation : 



A!i = hA-h-^ + hdzh-^ 
S^zw{g\ A) = S^zwig^; A^) - S^zwQih*- A^), with (41) 

g^ = hgh* 

for any element h of the complexified group Q'" . 

If we parametrize the impurity connexion as in previous section by A- = G^^djG with 
G G this action action has been cooked up such that : 

S^zM;A) = S^zA9'')-S^zAGG*) with A^ = G'^ ckG = - Al (42) 

This follows from eq.(l41|) and the fact that A^ = 0. 

This has two consequences. First we can explicitely evaluate the partition function for 
any sample of the disorder. Namely, we have : 

Z[A] = J Dge-'--^>^-'^^ = exp{S^zu,{GG*)) (43) 

In the last equation we used the fact that for the Haar measure on Q we have : / dgf{gigg2) = 
J dgf{g) for any gi, g2 G Q'^ and any function regular enough. 

The second consequence concerns the correlation functions of the ^-scalar local operators 
at fixed disorder. For example consider the operator 

i 

where (p^^iz) the primary fields in the representation R. These operators are invariant under 
the global Q symmetry. In the path integral formulation, they are represented by the insertion 
of tr{pji{g)), ie : 

<^R{z,z) = tr{pn{g)) (44) 
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where pr denotes the representation R oi Q. We then have 



{tr{pRig) ■ ■ ■)a 



Dg 

Z[A] 
Z[A] 



where we have used eq.(^2|). Now we see from eg. (^31) that the partition function Z[A] cancels 
against exp{Swzw{GG*)) . Therefore, using the invariance of the Haar measure, we obtain : 

MpRig)) ...)a = j Dg e-'--^^^ tr{pR{g)pR{GG*)-^) ■ ■ ■ (45) 

We have used the cychcity of the trace in order to reconstitute the product GG* . In other 
words, we have used the fact that the operator tr{pR{g)) is Q invariant. This is important 
because it proves that the G-dependence of these correlation functions projects from Q*-^ to 
Q'^/Q. Indeed, let us factorized G G as : 

G = HU with Heg^/g and Ueg (46) 

Then since UU* = 1, we have : 

MpRig) ■■■)a = j Dg e-^--^3hr{pR{g)pR{HH*)-') " " • (47) 

We thus have proved that the correlation functions of the neutral local operators factorize into 
the WZW correlation functions times an explicitely known function of the g'" /g component 
of the impurity connexion. This is analogue to the factorization of the previous section. 

Using the factorization (EH) we can now average over A. As before, it is convenient to 
change variable from A to G = HU. The Jacobian is given in eq.(^). In particular at 
a = oo, we obtain : 

{tr{pR{g) ■■■)a--- = I DHDg e^h^^^-^-^^^*)-^-^-^^) tr{pR{g)pR{HH*)-') ■ ■ " (48) 
Therefore, we have proved the following equivalence, valid in the ^-neutral sector : 



WZW^^\ = {WZWy^ X (WZW^c~/g'''^'j (49) 

where J\f is the number of correlation functions we are averaging. These are conformal field 
theories. The conformal dimensions of the operators in the disordered theory can be evaluate 
using the formula (pl|). 



4 The supersymmetric method : the random bond 
Ising model. 

In this section, we apply the supersymmetric method to study the random bond Ising model 
and the closely related Dirac theory coupled to a random potential and a random mass. 
These models have been analysed using the replica method in ref. ^ 
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As we will soon explain, the randomness of the bond interaction has a marginal effect on 
the large distance behavior of the Ising model. This allows us to use renormalization group 
techniques to study the infrared behavior of this random model. At criticality, the disorder 
only induces logarithmic corrections to the pure system. For example, the averages of the 
spin correlation functions at criticality are Q : 

N/A 



(a(i?)a(0))^~(^-j (log(i?/a))^ 
These have to be compared to the pure correlation functin which are : {a{R)a{0))^ ~ 

(*)"" 

Close to criticality, the specific heat of the Ising model possesses a logarithmic divergency : 
Cy = log(r) with T ~ -^^^H. In the disordered model, the behavior of the specific heat near 
criticality reads : 

Cy ~ — log (l + galog{-] 
ga \ r 

where ga represents the strengh of the disorder. This log(log r) behavior has to be compared 
to the (log r) behavior in the pure system. 

The perturbative study of the Dirac theory with a random potential and a random mass 
is very similar to the perturbative study of the random mass Ising model. However, the 
crucial difference is that contrary to the randomness of the mass, the randomness of the 
potential is marginally relevent. 

Although the supersymmetric method is only applicable to models which are gaussian 
at fixed disorder, it is certainly a good starting point for disentangling properties of a large 
class of disordered conformal field theories. The study of random bond Ising model clearly 
reveals that an appropriate algebraic framework for studying gaussian disordered systems 
at criticality should be based on affine Lie superalgebra with zero superdimension. In the 
case of the random Ising model, this algebra is OSp{2N\2N). The fact that the algebra has 
zero superdimension, ie. has an equal number of bosonic and fermionic generators, ensures 
that the Virasoro central charge of the Sugawara stress-tensor is zero. This is needed by 
construction for a disordered system. The vanishing of the central charge does not imply 
the triviality of the theory since it is non-unitary. The general framework should probably 
be based on the Wess-Zumino-Witten models, or their cosets, on Lie superalgebras which, 
like OSp{2N\2N), have equal numbers of bosonic and fermionic generators. The relevence 



of affine Lie superalgebras was independently realized in ref.pl 



•) The model. 

In its scaling limit (near criticality) the Ising model is described by a massive real Ma- 
jorana fermion with mass m ~ r = '■'^^"'^■^ with Tc the critical temperature. In presence 
of disorder, the mass becomes a function of the space position. The random Ising model is 
defined by the action {z = x + iy) : 

S[m{x)] = J — — (i>dzil> + ipdzip + im{x)ipip^ (50) 

where if) are grassmanian variables. The mass m{x) is coupled to the energy operator : 
e{x) = iip{x)ip{x). It is choosen to be a random quenched variable with a gaussian measure : 



P[m] = exp [ [ (m(x) — m)^ ] 

\ Ag J 27r J 



(51) 
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The energy operator e{x) has dimension one. The Harris criteria thus tell us that randomness 
in the bond interaction is marginal in the 2d Ising model. As we will see, it is not exactly 
marginal but only marginally irrelevent. 

Since, at fixed disorder the action (0) is a free quadratic action, we can use the su- 
persymmetric formulation to study the quenched model. In order to compute averages of 
product of correlation functions, we need to introduce an arbitrary number of copies of the 
fermions and of their supersymmetric partners. For simplicity, let us restrict to the case of 
two copies. The generalization to an arbitrary number of copies will be given later. Let us 
denote by ipi and ip2 the two real Majorana fermions, and introduce the complex fermions 
V'i by: 

^± = ^(^i±#2) (52) 

At fixed disorder, the fermionic action is: 

Sf = Si[m{x)] + S2[m{x)] 

^ + + ^^(V^-^+ - ^-?+)) (53) 

By definition this action can be used to compute products of two Ising correlation functions. 
Namely : 



{A{z) ■ ■ ■)m{x){B{w) ■ ■ = \ / Di)iDi)2 Ai{z) ■ ■ ■ B2{w 

Z\m[x)\ J 



where (■ ■ ■)m{x) refers to the Ising correlation function at fixed desorder. Ai{z) and B2{w) 
are the expressions of the Ising operators A{z)^ B{w) in the first and second copies. The 
partition function is a determinant : 



Z[m{x)] = I e-^^ = Det ^ 2 J ^ = n^rac (54) 

Its inverse can be represented as a path integral over bosonic complex fields t] and 7 : 

n{x)\ J 



Z[m{x)] 



with 



Sb = J ^ (^V^Tf + vdzl + i^^^^{V'y -V^)^ (55) 

Since the Dirac Hamiltonian H^,, as defined in (^l)), is purely imaginary. Hi = -Hd, the 
integral over the bosonic variables is an integral of an imaginary gaussian if we defined the 
complex conjugasion by r/* = 7 , 7* = rj. To insure absolute convergence, we could add the 
term e / ^(//7 + ?77) in the action 5*^. 

The total action is Stot = Sf+Sb- By construction, products of Ising correlation functions 
at fixed disorder can be rewritten as : 



(A(2;)---)„(,)(S(w) ■■■)„(,) = J D^±Dr]D-f A,{z) ■ ■ ■ B2{w) 



g Stot 
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The action Stot is supersymmetric. The supersymmetric transformation can be defined 

as : 



Si(j_ = 7] , Sipj^ = 
61] = , 57 = —■?/'+ 

Similarly for the ■ ■ ■. It is easy to check that the kinetic term in Stot is susy invariant: 
S{ip_dzip+ + V^^l) = 0- Moreover the field coupled to the mass terms in Stot is susy exact: 



This insures that the partition function defined by Stot is independent of the random mass 
m{x), exactly as BRS symmetry ensures gauge invariance in gauge theories. We actually 
have a stronger result : Stot is susy-exact. This is similar to topological theory; however, 
the difference with topological theory is that we are interested in correlation functions of 
non susy-closed operators. As we will explain below, this choice of the supersymmetric 
transformation is not unique. 

•) The effective action. 

We now set m = 0. I.e. we consider the effect of the disorder at criticality. Averaging 
over the quenched variables with the gaussian measure (|5l|) leads to the effective action : 



/d'^x / \ Q f d?'x 



S. + ll— %ert (56) 



with 



^pert = {lp-i'+ - i'-'ip+ +V1- Vl) (57) 

= 2('0]^?/)i)(?/)2'^2) + {susy partners). 

Clearly this action is supersymmetric since the total action Stot was supersymmetric for any 
values of the disorder before averaging. 

This action can be viewed as a perturbation of the (non-unitary) conformal field theory 
specified by the action 5**. This fixes the normalization of the fields to be : 

1 1 

, j{z)r]{w) ~ (58) 

z — w z — w 

The Virasoro algebra is the standard one. The Virasoro central charge is zero. Note that 
since the fermions tp± have dimension half, the perturbing field $pert has dimension two. It 
is therefore marginal. 

The conformal field theory specified by S"* is invariant under a supersymmetric algebra 
whose conserved currents are : 

G±(z) = r/(^)^±(^) , Gi{z) = 7(^)^±(^) 
K{z)=v\z) , K{z)=^\z) 
J{z) =: ij.{z)M^) : , H{z)=:j{z)r]{z): (59) 
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There are four fermionic currents, which are generators of supersymmetric transformations, 
and four bosonic ones. They form a representation of the affine 05*^(212) current algebra at 
level one. The root system of OSp{2\2) is given below : 



"^4 



K 



The black dots refer to the bosonic generators, while the symbols ® represent the fermionic 
roots. The non-trivial operator product expansions of the currents are the following : 



J{z)J{w) 



1 



J{z)G±{w) 
H{z)G±{w) 
H(z)K(w) ~ 



±1 



{z — w) 

1 

{z — w) 
2 



{z — wY 



G±{w) 
-K{w) 



{z — w) 

G±iz)G^iw) 



H(z)H(w) ~ 



-1 



{z — wY 

J{z)G^{w) ~ -^1-G±H 
[z — w) 

H{z)G±{w) ~ , ~^ M ±{w) 



H{z)K{w) 



+ 



{z — w) 
-2 



(60) 



{z — w) 



K{w) 



1 



G-{z)G+{w) 
K{z)G±{w) - 



1 



[z — w 
-2 



K{z)K{w) 
-K{w) 



(z — wY (z — w) 
2 4 



{H{w) ± J{w)) 



{z — w) 



G±iw) 



{z — wY {z — w) 
G-{z)G^H 

K{z)G±{w) ' 



H{w) 



{z — w) 
2 



K{w) 



{z — w) 



G±{w) 



The perturbing field can be written in terms of these fields : 

^pert = 2 



JJ-HH + ^{KK + KK) + - G^G+ + G+G^ - G+G^ 



In other words, the perturbation ( |^ is a current- current perturbation. It preserves the 
05*^(2 1 2) symmetry. 

The generalization to 2N copies is obvious. We introduce 2N copies of Majorana fermions 
"00, a = 1, • • • , 2N , or N copies of Dirac fermions ^Z'^, a = 1, ■ ■ ■ , A^, and 2N copies of their 
supersymmetric partners. The perturbing field is similarly obtained by a gaussian integral 
from the free random theory : 



^pert = 2 XI i'^ai^a){ipi3i^i3) + {susy partners) 



(61) 



The symmetry algebra is now extended to the affine OS'p(2A^|2A^) current algebra. The 
perturbing field is still bilinear in the currents. Hence, the random Ising model is described by 
a Gross-Neveu model on the Lie superalgebra OSp{2N\2N). In particular it is an integrable 
model, with factorizable scattering, Yangian symmetry, etc ... 
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•) Re normalization group computation. 

The action S^ff describes an interacting tlieory. Altliougli tlie perturbing field has di- 
mension two, and therefore is naively marginal, it breaks conformal invariance. To see it 
requires computing the beta function. As explained in the Appendix, the beta function at 
one loop is encoded in the operator product expansion (OPE) of the perturbing field. Since 
5* is a free gaussian theory, this OPE is easily computed using Wick's theorem and eqs. 
We obtain : 

g 

^pert{z)^pert{w) = 1 ^pert{w) + irrelcvent terms. 

\z — 

Using the formula ( |116D of the Appendix, we get the beta function : 

9 = Pig) = -9' + --- (62) 

The dots refer to higher loop corrections. We recognize the beta function for asymptotically 
free theory. It is easily integrated, giving the coupling constant flow : 

ga 

9r 



1 + ga \og{R/a) 

where Qa is the value of the coupling constant at the lattice cut-off a. Since Qa is positive 
by definition (it is the strengh of the disorder), qr decreases at large distances R. In other 
words, the theory is asymptotically free in the infrared regime. 

The beta function is independent of the number of copies, ie. of A^. 

•) Quenched correlation functions. 

The fact that the random Ising model is asymptotically free in the infrared regime is 
important since it allows us to compute the large distance behavior of two point quenched 
correlation functions using the renormalization group, (using formula (|117|) of the Appendix 
for the anomalous dimensions). 

We recall the OPE between the energy operator and the spin field. The fusion rules are 
e X 0" = ex. More precisely, we have : 

— C 

e(z)a(w) = i(ijjijj)(z)a(w) = - — ^-^^a(w) + ■ ■ • 

\z — w\ 

with Cl^ = 1/4. 

Let us consider the quenched average of products of two point functions of the spin field, 
Ie. {a{R)a{0))^ . Consider first the case = 1. In the supersymmetric effective theory, 
the quenched correlation function {a{R)a{0)) is represented by the two point function of the 
spin field cxi in the first copy. Its OPE with the perturbing field is : 

%ertiz)ai{w) = ^0{W) + ■■■ 

\z — w\ 

where 0{w), an operator of dimension (1/8 + 1), is irrelevent compared to ai which has 
dimension 1/8. Therefore, the anomalous dimension of cti at one loop is 7*^^-' = 1/8 + ■ ■ ■, up 
to irrelevent terms. Hence: 



1 \ 



1/4 



{a{R)am ~ ( -j (63) 
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It has the same infrared behavior as the spin-spin correlation function in the pure system at 
the critical temperature. 

Consider now the general case arbitrary. In the supersymmetric effective action, the 
quenched correlation function (cr(i?)cr(0))^ is represented by the two point function of the 
operator 0^'^\z), product of the spin fields in the N first copies : 

When computing the OPE between the perturbing field $pert and 0^'^\z) only the first term 
in eq.(|6l|), 2J2a<i3{4'a4'a){4'i34'i3), gives a relevent contribution. When computing the OPE 
between this term and 0^^\z) we can contract any pair of products of {il)^tl)a){i^ pi^p) with 
the pair of product of spin fields a^ais. E.g. 

(^iV'i)(^2V'2) ai{z)a2{z) ■ ■ - aN^z) 

Thus we have : 

^„ert{z)0^^\w) = _ ^(^ - ^'>^^'^- o^^\w) + irrdevent terms. 

\z ~ w\'^ 

Therefore, the anomalous dimension of 0^^\w) is 

-im N N{N-l)Cl 
' 8 4 ^ 

Recall that C^^^ = 1/4. It implies : 



(a(i?)a(0))^^^ = exp (^-2 J^^ ^da ) (a(a)a(O))^ 

N/A 



9R 



- (I) (log(/?/a))-(--)/« (64) 

There are logarithmic corrections to the pure system. 
•) The specific heat near criticality. 

In order to compute the behavior of the mean specific heat near the critical point we 
consider a non zero value of m: 

_ Tc-T 

m = — — = r 

Tc 

The specific heat Cy per unity of volume C„ = (If) "^^^ expressed in terms of the 
connected correlation function of the energy operator e{x) : 



a = / d'x {e{x)e{Q)Y 



Near criticality, its behavior is given by the infrared singular behavior of (e(x)e(0))^°"'"- at 
zero mass. The infrared cut-off is specified by the mass scale 1/m. Thus : 



l/m 

~ / d'x (e(x)6(0))™«" (65) 
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In the supersymmetric effective theory, the connected two-point correlation function for 
the energy operator is represented by the two-point function of the operator e'{x) = ei(x) — 
€2(3;) where 61^2(3^) refer to the energy operators in the first and second copies : 



2(e(x)e(0))~ = (e'(a;)e'(0)) 

To evaluate the behavior of C^, as m — requires evaluating the infrared behavior of 
(e'(x)e'(O)). This can be done using the renormalization group by computing the anomalous 
dimension of e'{x). The OPE between e'{x) and the perturbing field $pert follows from the 
fusion rule, e x e = 1 + e, in the Ising model : 

2 

<!>pert{z)e{w) = -e'{w) + ■■■ 

\z — 

Formula ( |117| ) of the Appendix gives : 



a^2 



Therefore, at large distances we have 
This gives the critical behavior of the specific heat : 

l/r 



a ~ / d'x{e{x)e{0)y 



~ -\og(l + gJogi-)) (66) 

This log(logr) behavior has to be compared to the (logr) behavior in the pure system. 
•) Miscellaneous remarks on the random Dirac theory 0. 

A model very closely connected to the random bond Ising model has been introduced in 
connection with the quantum Hall transition [0. It is a model of Dirac fermions coupled 
to a random potential, or more generally, to a random vector potential, a random mass and 
a random scalar potential. Its action is : 
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S = j ^{^lj-{d^ + A^)^+ + ^_{d, + A,)^^ 

+^^(^„^+ - ^_^_,) + + ^.tfj^)} 

The random variables A, m and V have a gaussian distribution with width qa, Qm and gv- 
Let us denote by $a/ and $y the perturbing fields coupled to the constants qa, Qm and 
gv after averaging over the disorder. In the two copy sector, we can write them in terms of 
the 05*^(212) currents as follows : 

^A = (H-J){H-J) 

^{<^v + '^m) = 2G^G+-2G^G^ + KK + KK 

i($y-<l>Af) = 2HH -2JJ + 2G^G+-2G+G^ 



These remarks arised from discussions we had with Martin Zirnbauer. 
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It is an interesting exercise to compute the beta function at one-loop for these couphng 
constants. We find : 



(3a = 9a = 32 gu gv^ 

f^v - Pm = gv - gu = 8 {gv + guf H (67) 

f^v + 13m = gv + gu = 8{gv + gn) (gv - gn) + 8gAigv + gn) ^ 

For generic value of the initial coupling constants this describes complicated flows. In the 
particular case with gA = gv = 0, we recover the beta function of the random bond Ising 
model. 

For gA = gu = 0, we have gy = 8gy. This means that a random potential is marginally 



relevent in the Dirac theory |19]. The coupling constant gy grows at large distances. The 
infrared fixed is expected to describe the quantum Hall effect transition. However, being a 
strong coupling problem no concrete description of it has been proposed. 

Note that in this perturbative computation we put all the epsilons which distinguished 
between advanced and retarded sectors equal to zero. In other words, we did not distin- 
guished whether we are computing averages of product of only advanced Green functions, or 
averages of product of advanced and retarded Green functions. It is likely that the infared 
behavior will not be the same in these two cases. It is tempting to conjecture that it will be 
trivial when we consider only advanced Green functions but non-trivial when mixed products 
of advanced and retarded Green functions are considered. 

Being formulated as current-current perturbation of a first order free theory, we can 
implement a Hubbard-Stratonovitch transformation on the random Dirac theory in order to 
produce a dual theory whose coupling constant is 1/g. The transformation goes as follows : 
one first represents the current- current interaction by a gaussian integral over an auxiliary 
gauge connexion A, and then integrate over the free fermions and bosons. This produces a 
dual theory whose fundamental fields are the connexion A, or equivalently the group element 
G such that A- = G^^dzG, and whose fundamental coupling constant is 1/g. However, in the 
mixed sector with both advanced and retarded Green functions, this naive transformation 
seems to be spoiled by divergences in the gaussian integrals. 



5 The replica method : the random phase sine-Gordon 
model. 

In this section, we apply the replica method to the random phase sine-Gordon model. It is a 
random version of the 2d XY model. It has also been used to describe other random physical 
systems; e.g. the 2d XY model in a random field p|, 0, interfacial roughening transition 0, 
randomly pinned flux lines in supraconductors P], etc... 

Let us first introduce the model from the XY point of view. Recall that the partition 
function of the XY model is defined as : 

It describes the thermodynamics of a fluctuating surface. The variables hi are interpreted as 
the height of the surface above a base plane. In the pure system the hi take integer values. 
The disordered model is defined by the same partition function but the height variables hi 
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take the values hi = di + rii with rii integers and di some random variables. Similarly as for 
the standard XY model, we can formulate this problem in an alternative form by using the 
Poisson summation formula : Zln/(^) = Sm/ (i$/(<l>)e*^'^"**. Introducing the variables 
one gets : 

Z^jD^^exp ( -| - ^jf + i2nJ2mj(^j - dj) 

In the continuum limit, the variables $j are replaced by a field $(a:) and the discrete Lapla- 
cian becomes a continuous one. Near the critical point, only the first harmonic of the local 
potential is relevent. Keeping only this first harmonics gives the action : 

/^(f(5.*)'-Acos($(a;)-d(a;))) 

with d{x) a random quenched fields. This is the model we will use to illustrate the replica 
method. 

As we explain below, this model as two different phase : a high temperature phase in 
which the disorder is irrelevent and a low temperature phase in which it is relevent. In the 
high temperature phase nothing interesting happens. In the low temperature phase, the 

large distance behavior of the system could a priori be different than in the pure system. 
However, a precise description of this behavior is still missing. For example, the behavior 
of the quenched correlation functions of $(a;) is still controversial. One finds that at large 
distance we have : 



( [^{x) - $(0)]2 ) = yl(log \x\) + BiT - r,)^(log \x\y 

However, renormalization group computations predict a non-vanihsing S 7^ 0, while varia- 
tional approaches give B — In the sequel we will present both approaches. 

•) The model. 

Let us first rewrite the action with more appropriate notations. We also need to slightly 
generalize it by introducting a random potential in addition to the random phase. The action 
reads : 

= 1^ (f (^'^^)' - ^'^^'^^ - ^(^)^'* - r(^)e-^*) (68) 

The constant K is proportional to the inverse temperature, K oc l/T. In the formula, A^{x) 
and ^{x) are random quenched variables with gaussian measure : 



P[A] = exp 
P[e] = exp 



1 f d'^x 
Yg 1 

1 f d^x 

/ 

2a J Att 



(69) 



In absence of disorder, the dimension of the vertex operator exp(i$) is : dim(e'*) = 
The Harris criteria then tell us that there are two different phases : 

low temperature phase {K > Kc — 1) =^ disorder is relevent 
high temperature phase {K < Kc — 1) ^ disorder is irrelevent 
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Recall that K ocl/T. At high temperature, the large distance behavior is identical to that 
of the pure system. At the critical temperature K = Kc, the disorder has a marginal effect. 
At low temperature, the large distance behavior is a priori different. 

•) Symmetries. 

Similarly as for the random gaussian model, the random phase sine-Gordon model pos- 
sesses a U{1) symmetry whose Noether current corresponds to insertion of 9^$ in the 
quenched connected correlation functions. This symmetry in particular implies that all 
the quenched connected correlation functions of (9^$ are unaffected by the disorder. 

Just as for the Gaussian model, we can decompose as = d^A + e^^d^(. Then the 
field ( decouples from the action and from the measure. So we can set it to zero, and we 
have : 



P[A] = exp 



2g J 4n 



id.AY 



with A^ = d^A (70) 



We will denote by S'($|A; ^ the action (|68|) with A^ = d^A. 

The simplest way to visualize this symmetry consists in introducing sources for 9^$ in 
the action (|68D : 

72 

5(<I>|A;0-5(<I>|J^;A;0 = 5($|A;0- / 

Let us decompose as = d^p + e^i^d^r]. The field r] decouples and only the field p is 
relevent. The U{1) symmetry follows form the fact that the p and A dependence can be 
absorbed into a shift of $. Namely : 

5($|J,; A,0 = ^("^ - ^I^M = ^ = O'O - ^/ ^{d,{p + A)f (71) 

with ^ = ^e*^''^^)/^ and = d^p. As a consequence, the generating function of the con- 
nected correlation functions of 9^$, which is log Z[J^, A, ,^], can be written as : 

logZ[J„A,e] = '^{d,{p + A)f + log Z[J, = A = 0,e] (72) 

Now, notice that ^ and ^ have identical distribution. Therefore, averaging ([7^ ) over the 
disorder implies : 



logZ[J^, A,e] = j ^(d^pf + const. 



where the constant is independent of J^. This shows that the quenched average of the 
connected correlation functions of the U{1) current are identical to the correlation functions 
of this current in absence of disorder. In other words, the connected correlation functions 
are protected from the disorder by the U{1) symmetry. 

It is worth noticing that the presence of the disorder restore the U{1) symmetry which 
is absent in the pure sine-Gordon model. 

In connection with this U{1) symmetry, the model possesses the remarkable property 
that the (yf-dependence of the correlation functions of the vertex operators can be factorized. 
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Consider eq.([7T|) in absence of source : = 0. The fact that the A-dependence can be 
absorbed in a translation of $ imphes for the correlation functions that : 

Let Gaj^^...{xi, ■ ■ ■ \g, a) be the quenched correlations of the vertex operator exp{ia^) : 



G^,,...{xu ■■ -19,(7) = (e-i*{-i) ■ ■ ■)A,g(- • (73) 

Integrating over A using the free field gaussian measure (|7D[), and using the fact that ^ and 
^ have the same measure, we deduce : 

...\g,a) = l[ \xp - ...\g = 0,a) (74) 

p<q 

Equivalently, 

dgGa^^-{xi,- ■ ■\9,0r) = (^'^^^^Og{\Xp- Xg\^)^ G„,,...(xi, ■ ■ ■ jfii, cr) (75) 

This identity will be useful for analyzing the renormalization group equations. 
•) The effective action. 

Since the action (|68D is not free we cannot rely on the supersymmetric method but only 
on the replica trick. Therefore, as explained in the introduction, we consider n copies of the 
system with the same disorder and then average over the disorder. This gives the following 
effective action : 



Let us first look at the kinetic term. It is of the form / ^^{dfj,^'')Grs{dfj,^'') with 

Grs = KSrs - g = {K ~ g)6rs - g{l - Srs) 

Since the interaction only involves the difference of the replicated fields, it is convenient to 
decompose the kinetic terms as : 



Note that the field {J2r ) decouples. Its correlation functions are therefore unaffected 
by the disorder. Since these correlation functions represent the averages of the connected 
correlation functions of the U{1) current, we recover the previous result obtained from the 
U(l) Ward identities. In particular, (9^(1],.$^) is a U{1) conserved current. We could also 
use this decomposition to derive eq. ([75|) . 

The kinetic term ( [77| ) fixes the normalization of the vertex operators. In particular, we 
find the dimension of the perturbing field : dim ^^e*'-*' = It is independent of a. It 
is relevent for > 1 : we thus recover the Harris criteria. 
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This decomposition has a Lie algebraic interpretation. It corresponds to the decomposi- 
tion of U{n) = SU{n) x U{1). The second term in (|77D can be rewritten as ^{d^i^Y where 
takes values in the Cartan subalgebra of SU{n). In terms of (p, the perturbing field in (|76D 
reads : X)a6xp(m.(^) where the sum extends over all the SU{n) roots. At K = K^, the fields 
exp{ia.ip) have dimension two, and the perturbing field is a current- current interaction. The 
action (|76D at K = Kc is therefore equivalent to the SU{n) Gross-Neveu model. We will 
later see that at this point the model is asymptotically free in the infrared. Hence at the 
critical temperature K = Kc, the disorder only induces logarithmic corrections. 

•) Re normalization group. 

The renormalization group allows to perturbatively analyse the behavior of the system 
in the low temperature phase. We will do a one loop computation, which is valid close to 
the critical temperature, ie. ^J^^" <^ 1- 

As we are already familiar with, the one-loop beta functions are encoded in the OPE of 
the fields. Let us introduce the following notation : 

O, = ^exp(z($^'-$^)) (78) 

O2 = (79) 

The field Oi is the perturbing field associated to the coupling constant a. The field O2 is 
one of the kinetic field. Notice that O2 = {idz0){id^0). It is necessary to introduce it since 
it is generated from Oi by OPE. Indeed we have : 

2(n — 2) 2n 
Oi{z)Oi{0) = -^^-^jj^Oi{0) + j-^^jj^O2{0)+irrelevent terms, (80) 

2/^2 

O2{z)Oi{0) = -^-n^OiiO) +irrelevent terms, 



\z 



O2{z)O2{0) = irrelevent terms. 

For n = 1 these are the familiar OPE of the Kosterlitz-Thouless transition. Note the change 
of behavior between n < 2 or n > 2. The first OPE indicates how cr is renormalized, but the 
second equation shows that the kinetic term needs also to be renormalized. 



Using the formula (|116|) of the Appendix, we get the beta functions 



Pa = 2 j^^" ) a + (n - 2) cT^ + • • • 
13k = 2^ + ■■■ 

This shows that the coupling to would have been generated at one loop even if we did 
not start with it. Since the field (X]r decouples, the coupling [K — ng) is unrenormalized 
(at any order in perturbation theory). Thus : Pk = n^g. Setting = as required by the 



replica trick, we get : 



Pa = 2(^^)a- 2a' + ..- (81) 

P. = y + 
Pk = Q 
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So, K is unrenormalized at ra = 0. It may appears surprising that the couphng g is renor- 
mahzed although we know exactly the (^-dependence of any correlation functions. We will 
later see that this is not in contradiction with the renormalization group equations. 

From equation (|8TD, we immediatly see that in the low temperature phase {K > Kc), the 
beta function P„ possesses a non trivial infrared zero at a* : 



K-K, 



+ 



for 



< 1. 



^2) 



Notice the fact the other beta function (3g does not vanish at cr* : 



I3g{(y*) 



1 fK-K, 



2\ Kr 



+ 



(83) 



Hence, even at a* the coupling g will continue to flow. We may characterize such pseudo- 
fixed point as a "run away fixed point" . This is a particularity of the model which has direct 
consequences on the correlation functions. 

•) Quenched correlation functions. 

Let us now analyse the renormalization group equations. As we already said, the field 
would have been generated at one loop if not present initially. But on other hand we know 
that the dependence on g can be completely disentangled. This apparent conflict should 
have an effect on the renormalization group equations. 

Consider the correlation functions Gaj^^...{xi, ■ ■ - {g, o") defined in eq.(^). They satisfy the 
renormalization group equations : 



da 



G'ai,...(a;i, ■ ■ ■ \g,(T) 







where jp{g,a) are the anomalous dimensions. But the (^-dependence is explicitly known, 
therefore using eq. (|75D we get : 



E 



x: 



d_ 



d 



+ ^ 7p(^7, ^) - f3<r{a)— - P,{a) ^ ^ log(|a;, - 



p<q 



Note that it is now possible to set g equal to zero. In particular at the infrared fixed point 
a^, in which /?o-(cr*) = 0, we get : 



E ^.^^ + E 7.(-*) - E ^ iog(k. 



G*a,,..{xi,---\g,cr) 



with (3g{(J^,) given in eq.(^3|). The effect of the g-flow is to add the extra logarithmic term in 
the renormalization group equations (|8^). 

This can be used to compute two-point functions at the infrared fixed point. Consider, 



G,{x) = (exp(^a($(x)-$(0)))^^ 
G2{x) = {exp{ia^{x))){exp{-ia^{0)))^ 



30 



Let 71,2(0"*) be their anomalous dimensions at the infrared fixed point. The renormahzation 
group equation (^4D gives : 

G,,2ix) = |a:|-2^^'^(-) exp (-^^^^(log |a:|)^) (85) 

Notice the (log|x|)^ correction which arises from the renormahzation of g. This term is 
independent of the anomalous dimension. 

The anomalous dimensions are even in a and vanishes at a = 0, therefore : 

71,2(0-*) = — Pi,2(o-*) + C(a^), with pi(cr*) = 1 + 0' 



Expanding in power of a^, gives the two point functions of $ : 

W^cT^WW) = log \x\ + §^(log \x\r (86) 



iHx) - m) r = log ixi + §^(iog \x\r (8?) 



Note that the (log |x|)^ cancels in the connected correlation function {[^{x) — $(0)]2)conn as 
it should be, since this connected correlation function is unaffected by the disorder. 

There is a crossover from a (log |x|) to a (log |x|)^ behavior. For |x| < Rcross "we have a 
usual (log|x|) behavior, while for |x| > Rcross we have a (log|x|)^ behavior. The crossover 

length is Rcross with 

log Rcross r,/ . ~ tr^ V^2 ' /^'^ {K-Kc)<^l. 

Rcross is exponentially large close to the phase transition. However eqs.(^,0) are true to 
all order in perturbation theory if we can rely on the renormahzation group in the replica 
symmetric approach. 

We can use eqs.(|^,^) to find an estimate of the width of the surface for a system of 
finite length L, ie. to find an estimate of ($(0)^)^. We define it as the integral of the Fourier 
transform T{q) of (|86D using 1/L as ultraviolet cutoff. At short momenta q -C 1/ Rcross, the 
propagator r(g) is dominated by the Fourier transform of the (log|x|)^ term; i.e. T{q) ~ 

/?.(^*)(^)- Thus: 



QiL) = {m')L = [ d'qT{q) ~ /3,(a,) (logL)^ 

Jl/L 



n/L 

for L ^ Rcross- This has to be compare with the pure case which gives a (-^ logL^) behavior. 
•) A large N model. 

The formula (|S^, |S^) are still controversial, theoretically as well as numerically. Vari- 
ational approaches, (part of which we will describe below), predict a (log|x|) behavior 
21\, B3, |2^|. The RG flows was also found to be unstable again assymmetric replica pertur- 



bations |3T[]. The numerical veriflcations of (^) are also not settled : a (log |a;|) behavior 
was found in ref. |^, while more recent simulations seem to indicate a (log|a;|)^ 
behavior. In view of this conflict, and since the variational approaches are argued to be 
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exact for systems with a large number of components |35|, in ref.[0 we studied a large 
version of the model (|68D. Using RG computations based on the (a priori symmetric) replica 
trick, we find that our large N model possesses a non-trivial infrared fixed point in which 
the correlation functions have the form [87| ) but with the (log |x|)^ term suppressed by a 
factor (1/A^'^) compared to the (log|x|) term. 

To introduce the large version of (|68D, it is convenient to fermionize it. The fermionic 
form of the random phase sine-Gordon model is a massless Thirring model coupled to a 
quenched potential and a random phase ^. To define its large N version, we need to 

introduce Dirac fermions with components ipj_ and ipj_ with k = 1, ■ ■ ■ , N. Let z = x + iy 
and z = X — iy he the complex coordinates on the plane. The action is : 

j2 / N N N \ 

^^^^ = /-^ E(V^-;^^++^-;fc^^V^+)-^(E^-;^V't)(E^^^^^ (89) 

^ \A:=1 k=l k=l I 

- / — mi: V'-;.^') + MH i^-A) + ^(E + r (E 

\ / 

In absence of disorder, it is conformally invariant. The random potential = {A-, A^) is 

coupled to the U{1) currents = J2k=i'4'-;k^+ and J- = J2k=i'^--k'^+ of the unperturbed 
theory. At ^ = 0, the random potential does not break conformal invariance. 

There are a priori many ways to generalize the action (^) to a large A^ version. The 
action (|89D has been designed in such way as (i) to keep the number of disordered variables 
fixed, (ii) to preserve the exact conformal invariance in absence of disorder, and (iii) to 
parallel as much as possible standard properties of large A^ models. 

The fermionic action (^) can be bosonized back using non-Abelian bosonization |37 . 
As usual, since the pure system describes A^ Dirac fermions, the pure bosonized theory will 
be described by a su{N) Wess-Zumino-Witten (WZW) model at level one plus a massless 
free field. The su{N) WZW model at level one possesses primary fields taking values in the 
(A^ — 1) fundamental representations of su{N). Let 0^ and (p^k be the chiral WZW primary 
fields which take values in the defining representation of su{N) and in its complex conjugate. 
Their conformal weights are both equal to ("^1^) • Let us denote by (p the gaussian free field. 
The original fermions can be written as the product of these WZW primary fields by a 
vertex operator of the gaussian model. Namely, ip'l = (p^ e^'^ and ip-.-^k = 'Pu,k e '^'^ , and 
similarly for the other chiral components ipj^ and '4'--k- The bosonic form of the action (pOD 
is : 

~ / 17 {M^^^uV + ^(^)'^ ""^^ + C{xW e"^^) (90) 

where S^zw refers to the su{N)i WZW action. We have introduced the composite fields 
^{z,z) = Ef=i0o(^)^(^) and $* = T.k=i(t^^k{z)^{z). Equivalently, $ = tr^{G) with G 
the group valued field of the WZW model. For A^ = 1 the WZW terms are absent and we 



recover the action (|68|) of the random phase sine-Gordon model. 

The renormalization group computation done in ref. |^ predicts a non-trivial infrared 
fixed point in the low temperature phase with correlation function of the form (R3) but with 



32 



(3g{(T^) given by : 

The occurence of this factor could explain why the (log|x|)^ term does not manifest 

itself in the variational approaches. 

6 The variational method : the random phase sine- 
Gordon model. 

Various variational approaches to disorder systems have been proposed. They can be applied 
after or before disorder averaging. In the first case, one first averages over the disorder using 
the replica trick, and then implements a variational method on the replicated model, cf. eg 
[ P^ . In the second case, one applies a variational method on any sample of fixed disorder, 
and then averages over the disorder. We will present the second method using the random 
phase sine-Gordon model as example, following ref.[^ . But it can clearly be applied to 
other models. The variational method at fixed disorder leads to an exact bound to the 
averaged free energy. 

Consider the partition function Z[,^(a;)] = / D$ exp(— 5'($|^)), where ^(^l,^) is the action 
( |68D at = 0, which we recall : 

Sim = / ^ (f (^^•^)' - ^(^y - r(^)e-^*) 

At fixed disorder, we approximate the partition function by a gaussian action : 

Some) = \j £xd% ^x)G-\x - vMy) = IJ^, mG-\k)m (9i) 

Here the hatted quantities refer to the Fourier transforms. Notice that the ansatz (^) is 
choosen to be translation invariant, and that the gaussian is centered around the origin. We 
choose a gaussian anstaz otherwise the computations are undoable. The kernel G{x — y) is 
the variational parameter. 

The approximated partition function is = / D$ exp(— 5*0). Let F and Fq be the 
respective free energies, ie. Z = and Zq = e~^°. For any realization of the disorder and 
for any choice of 5*0, the following inequality holds : 

F<Fo + {iS-So))o (92) 

where (■ ■ ■)o refers to the expectation values with the measure Sq. This is proved using the 
following inequality : 

Therefore, the best approximated action will be find by minimazing (Fq + {{S — So))o). 
That is : 

^ (Fo + ((S-5o))o) =0. (93) 



6G{k) 
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The inequality (p2^ gives a upper bound to the free energy. 

Let us apply this general setup to our example. As it is formulated the method is more 
appropriate to study the free energy than the correlation functions. Thus we consider the 
system in a box of finite volume of size Vol. = L"^ and look for the L-dependence of the free 
energy. That is, we use the variational free energy to analyse the finite size effects. 

Since 5'o is Gaussian, Fq and {{S — So))o are easily computed. The free energy Fq per 
unit of volume is given by the logarithm of a determinant : 



1 



d'k \ogG{k) 



Vol. 2 

while for the expectation value of {S — So) per unit of volume we have : 

{{S - 5o))o _ .K f d^k Q^^^ _ r fx^^^^^ ^ ^.^^^^ ^^^^ 



Vol. 2 J (27r)2 ^ ' Vol. J 47r 

with 

Q = ($^(0)) = G(0) = / (0G(A:) (95) 

This parameter is naturally interpreted as the width of the system. 

Notice that the expectation value (0) only depends on a particular moment of ^(x). It 
does not depend on all the details of the disorder configuration. It is this particular fact 
which makes the variational approach doable. 

The variational equations (|93D are simple to compute. They determine the kernel G{k) : 

G{k) = % f zr^) (96) 

where the effective mass is : 



K Kk'^ + M'^ 



We have introduced the prefactor (Vol.y^'^ in the definition of Dq to make it scale invariant. 
Since the effective mass depends on Q, which is a functional of the kernel G{k), eqs. (|95| , 



96) form a set of non linear coupled equations for G{k) which we rewrite below : 



The effective mass depends on the disorder through its moment Dq. This moment can 
be either positive or negative, with a symmetric probability distribution around the origin. 
Thus we have to study separately the two cases : Dq > or Dq < 0. This analysis was done 



in ref. 29 



For Dq > 0, the effective mass is real and the Green function G{k) has no pole. Therefore, 
for Q we get : 

/•V'^ d^k 1 + ]Vf2\ 
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where kuv = 1/a and kj^ = 1/L are the ultraviolet and infrared cutoff. Recall that the 
volume is Vol. = L?. Inserting this expression in the definition (|9^) of the effective mass 
gives : 

The effective mass vanishes as L ^ oo, but with different power of L for i^' < i^'c = 1 or 
K > Kc = \. These differences arise from the different behavior of the last term in eq. (|100D . 
For K < Kc, we have <C (1/L)^, and therefore the last term in eq.( |100| ) behaves like 
For K > Kc, the mass term dominates, ^ and therefore, the last 

term in ( [LOOP behaves like M^^^ . Hence, for Dq > we obtain : 

M\L) ~ Qy^" ; Q(L) ~ -^logL^ for K < (101) 

(M2(L))1-w ^ (^ij ; g(L)~^^L_iogL2 for K > (102) 

For Do < 0, since the situation is quite different since the effective mass square is negative. 
The Green function G{k) now possesses a pole at fc^ = — M^. In order to analyse the effect 
of this pole we have to remember that for a system in a box of length L, the momenta are 
quantized to discret values : (kx,ky) = (^^,^^), with {nx,ny) integers. So in the key 



equation (p8|) the integral is actually a discret sum : / j2J2n^,n^- In eq. 



we 



seperate the first terms which correspond to (n^. = ±, rij^ = ±) from the others which we 
approximate by an integral. We obtain : 

4:71 4 

Q{L) = ^-—2 + 

+ llog(LVa2) (103) 



Att 
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where we have neglected the dependence in the last integral. Once again, the effective 
mass behaves differently as L ^ cxd for K < Kc and K > K^. These different behaviors are 
distinguished by the relative importance between the two terms in eq.( |103D . Indeed, suppose 
that ^ 1/L^. Then the first term in ( |103|) is irrelevent and therefore Q{L) ~ -^logL^. 



However, inserting this value of Q in the definition of as a function of Q, cf eq.( 
we deduce then that M'^{L) ~ {1/LY~^k. Thus consistency of the hypothesis <C 1/L^ 
requires K < 1. When becomes of order 1/L^, the first term in ( 103 ) dominates and 



remains frozen to this values. The expression of as a function of Q then tell us that 
Q{L) ~ logL^. Summarizing, for Z^o < we get : 

M2(L) ~ (^iy^"" ; Q(L)~-ilogL2 for K < K^ (104) 

M2(L)~Q)' ; Q(L)~logL2 for K > K^ (105) 

Notice that in the high temperature phase K < K^., the behavior of Q{L) is the same for 
positive or negative, while this behavior is different in the low temperature phase. 
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These behaviors have been obtained at fixed disorder. We can now average over the dis- 
order. Since the effective mass was only a function of Dq which is symmetrically distributed 
around the origin, the average value of Q{L) is half of the sum of its values for Dq positive 
and negative. Hence, at large L, we have : 

far K <K^ (106) 
for K> (107) 

The averaged behavior at high temperature is the same as in absence of disorder. That is, the 
disorder is irrelevent in the high temperature phase as we found in the previous section using 
the replica approach. But the disorder is relevent in the low temperature phase. However 
the result obtained for Q{L) in the low temperature phase disagree with the result obtained 
sing the symmetric replica trick, cf eq.(^). 

This variational method gives poor results for the correlation functions. Indeed, since 
we choose the Gaussian anstaz ( [9T| ) to be centered around the orgin, the variational one- 
point functions vanishes : ($(x))o = 0. But the connected correlation are unaffected by 
the disorder since it is protected by the U{1) symmetry. Therefore, the variational two- 
point function is also unaffected by the disorder, which is probably not realistic. A more 
appropriate ansatz could be to choose a gaussian action not centered around the origin. 

Also, this variational approach does not take one-loop effect into account. This could 
be the origin of the disagrement between the renormalization group and the variational 
approaches. 



g(L)~-logL2 



Q(^) 



K 



2K - 1 



loffL^ 



7 Replica symmetry breaking or not? 

In this section we very shortly describe how replica symmetry breaking is incorporated in 
the renormalization group perturbative approach based on the replica method. The aim is 
not to present all the details and subtilities of the replica symmetry breaking, (there already 
exist extensive reviews on the subject), but only to introduce the main steps. 

As example we choose the minimal conformal models perturbed by random bond inter- 
action. The replica symmetric was studied in ||38| , while the study of the theory with replica 



symmetry breaking was done in |39]. The basic examples are the random bond Ising or Potts 
models. 

The simplest minimal model is the Ising model, whose disordered version was studied 
above in the supersymmetric approach. In the scaling limit, random bond interaction is 
represented by a perturbation by the energy operator e{x) with a random coupling constant. 
In the Kac classification this operator is the $12 operator. Thus we consider the random 
models : 



S[g{x)] = S* + J d^xg{x)<^i2{x) 



where S^, represents the action of the corresponding minimal conformal model, e.g the Ising or 
Potts models. As we have seen, in the Ising model the disorder is marginal and only induces 
logarithmic corrections. The dimension of $12 in the Potts model is dim($i2) = | < 1. It is 
therefore a strictly relevent disordered perturbation. 
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This models were studied in ref. using perturbative renormalization group cpmputa- 
tion in the a priori symmetric rephcated theory. 
After rephca, the effective action is : 

Seff = T.S: + ^ [ d'xY^^UxWAx) (108) 



In eq.( |108|) , one chooses to restrict the sum to r 7^ s since in the OPE of two $12(2:) operators, 



which is given by the fusion rule $12 x $12 = 1 + $13, only the identity arises with a singular 
term. The compatilibity of this hypothesis with the renormalization group has to be checked, 
and this is not a priori clear. 

Replica symmetry breaking is incorporated in two steps One first promotes the 



coupling constant cr to a matrix cr^s , and then consider ars at n = has a hierarchical Parisi 
matrix. This amounts to parametrize it by a diagonal element and a function (t{x) with 
X e [n, 1], n — *• : 

a ^ars^{?f,(T{x)) (109) 

The multiplication law of two Parisi matrices parametrized by (ai,(Ji(x)) and {d'2,<J2{x)) is 
then defined by : 



with 



(ai, aiix j) ■ (?2, (T2(a;)) = {h, h{x)) (110) 

h = aid-2-d\d^ (111) 
h{x) = -nai{x)a2{x) + {ai -ai)a2{x) + {a2 -W2)(J2{x) 

- r dy{ai{x) - ai(y))(a2(x) - a2{y) (112) 



where f = Jn dxf{x). 

One can now study the consequences of this anstaz in the renormalization group. This is 
done by first computing the renormalization group equations with the matrix ars and then 
implementing Parisi's ansatz. As usual the one loop beta functions (5rs = c^rs are encoded 
in the operator product expansions. However, non replica symmetric fixed points in the 
random mininal conformal models only appear at two loops. 

The two-loop beta function was computed in ref.||3^. using an epsilon expansion. The 
central charge is parametrized as : 

c = 13 — 6(a^ + alf;^) with = e 

3 

The case e = corresponds to the Ising model c = |. The Potts model corresponds to 
e = The dimension of the $12 primary field is then : 

3 

dim($i2) = 1-2^ 

For e = the disorder is marginal, while for e <^ 1 the disorder is slightly relevent. This 
allows to implement an e-expansion, and in particular to determine the non-trivial fixed 
point in an e-expansion. Since we assumed that the diagonal matrix elements of ars vanish. 
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in Parisi's ansatz the matrix ars is represented only be a function ct{x). The vanishing of the 
beta functions then determines the possible fixed point function a* (a;). According to ref . ||39|| , 
the fixed point equation reads : 



3ecr*(a;) - 2(j^,(j^{x) - / dy{(j^{x) - (T^{y)Y + crl{x) + cr^a^^x) = (113) 



with a* = Jq dxa^{x). 

The replica symmetric solution corresponds to cr*(a;) = const 



3 9 



a^:{x) = const. = -e + -e^ + ■ ■ • 
2 4 

The non replica symmetric solution is given by a solution for which a^{x) growth linearly 
for < a; < Xi and then remains constant for xi < a; < 1. This is called a one-step replica 



symmetry breaking. It is found by deriving eq. (|113|) with respect to x as many times as 
necessary. Explicitly, one has : 

, . f ^x, if < X < Xi, .,, 9 27 2 

ct^kIx) = < y -r ^ With Xi = -e H e + ■ ■ ■ . 

^ ' \ \xi, if xi < X < 1. 2 2 

This solution is not present at one-loop. It has been checked that it is a stable solution of 
eq-dig). 

The symmetric and non-symmetric solution can be distinguished by analysing the anoma- 
lous dimensions at the corresponding infrared fixed points. For example, the dimension of 
the energy operator at the infrared fixed point differs in the two solutions by two percent up 
to 0{e^). There is up to now no numerical evidence in favor of the non-symmetric solution 
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8 Appendix: Renormalization group and OPE. 

In this appendix we gather standard informations about the renormalization group in two 



dimension, cf eg. ||43[| . 



Consider partition functions and correlation functions computed with the measure / Dcf) exp( 
with a perturbed action : 

S = S, + J2g' [ d^x^iix) (114) 

where $i(x) are relevent primary operators of dimension hi. Suppose that these fields satisfy 
the following operator product expansion : 



Then, the beta function at one loop is : 



= (3\g) = (2 - K)g' - tt J] g' + • • • (116) 

jk 
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No summation in the first term but summation over j, k in the second. The summation is 
over all the relevent fields generated by the product operator expansion. 

In the same way, if Oa is a set of operators with OPE with the perturbing field given 
by the structure constant C^^,. Then the matrix of anomalous dimensions (7 = —adalogZ) 
is : 

i:,=hJ:, + 27:Y.Cf^'9' + --- (117) 
j 

Note that we have : 7] = 26} - d^pK 

The renormalization group equations are : 



{0,{x^)---Op{xp))^Q (118) 
For the two point functions, the integrated version of the RG equation reads : 

{0{R)Om,^a) = (O(a)O(0)),(«) exp (-2 f^^^^ ^^^) (^S) 
where 7$ is the 7-function for 
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